Optimizing entangling quantum gates for physical systems 
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Optimal control theory is a versatile tool that presents a route to significantly improving figures 
of merit for quantum information tasks. We combine it here with the geometric theory for local 
equivalence classes of two-qubit operations to derive an optimization algorithm that determines 
the best entangling two-qubit gate for a given physical setting. We demonstrate the power of this 
approach for trapped polar molecules and neutral atoms. 
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I. INTRODUCTION 

Quantum information science requires extreme accu- 
racy in implementation of quantum tasks such as gate 
operationsLin order to fulfill the stringent criteria of fault 
tolerance [Ij. However the dynamics underlying a gate 
operation is often very complex, occurring in a Hilbert 
space much larger than that of the qubits. Optimal quan- 
tum control then provides an indispensable tool for ob- 
taining a high-fidelity implementation Realization 
of a universal set of gates comprises multiple levels of 
difficulty, because different physics is typically involved 
for one- and two-qubit operations. For many qubit sys- 
tems, two-qubit gates are the most challenging since they 
involve controlled interaction between two otherwise iso- 
lated quantum systems. Furthermore, for a given phys- 
ical implementation, it is not necessarily a priori clear 
which two-qubit gate can best be implemented when the 
practical criteria of optimal achievable fidelity and real- 
istic gate operation time are imposed. 

We address this problem here by combining optimal 
control theory [3| with the geometric theory for local 
equivalence classes of two-qubit operations 3]. This al- 
lows us to develop an algorithm that optimizes for the 
non-local content of a two-qubit gate rather than for fi- 
delity of a specific gate such as CNOT. The resulting 
separation of non-local from local control objectives re- 
laxes the control constraints and enables both maximum 
fidelities to be reached and fundamental limits for control 
to be identified. We apply our algorithm to trapped po- 
lar molecules and neutral atoms, both candidates for re- 
alizing quantum computation. Manipulation of trapped 
polar molecules with microwave fields has been shown 
to allow realization of effective spin-spin models [3] with 
continuously tunable parameters. In our first example we 
use the combined optimal/geometric control algorithm 
to determine which non-local gates can be realized for a 
given underlying molecule-microwave field Hamiltonian. 
In our second example, we obtain optimal solutions for 



a Rydberg gate [a] between trapped neutral atoms in 
the presence of both decay processes and entangling cou- 
plings between internal and external degrees of freedom. 
We show that high fidelity Rydberg gate operation is pos- 
sible even for a configuration where the blockade regime 
is not reached and despite spontaneous emission from in- 
termediate states. 



II. OPTIMIZING THE NON-LOCAL CONTENT 
OF A TWO-QUBIT GATE 

A. Optimal control theory 

High-fidelity implementations of quantum gates can be 
obtained with optimal control theory by defining a suit- 
able distance measure between the desired unitary O and 
the actual evolution, e.g.. 



Ji? = l-l9^e 



Tr{6+PArU(T,0;e)P^} 



(1) 



and minimizing it with respect to some external field 
e{t) [Ij- Here, U(T, 0;e) represents the evolution of the 
system under the action of an external field from time 
to time T. For example, e{t) can be a pulsed laser field 
or a time-dependent magnetic field. The Hilbert space in 
which the system evolves is possibly very large. The log- 
ical subspace, i.e., the subspace of the total Hilbert space 
in which O acts has dimension N, N = d\mHo = 4 
for a two-qubit gate). Pjv denotes the projector onto this 
subspace. The trace is evaluated by choosing a suitable 
orthonormal basis of the subspace Hoi {\^k=i,...,N)}- 
The evolution of the system is thus expressed in terms of 
the time evolved basis states, \(pk{T)) — 0{T,0;e)\Lpk), 
and becomes a functional of the states \(pk{T)). For 
the specific choice of Eq. ([1]) , is a phase sensitive func- 
tional that depends linearly on the states, \tpk{T)), and 
equals zero for perfect implementations of U. 
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Additional constraints can be introduced to ensure fi- 
nite pulse fluence, 



9a = K [e{t) ~ ErefitT / S{t)dt , (2) 

Jo 

or avoid population of states subject to loss, 

9b^J^ ^{'Pmit)\Pavoid\'^niit))dt. (3) 

■^0 m=l 

Here, eref{t) denotes a reference field, S{t) is a shape 
function to switch the field smoothly on and off, and Xa 



and Ab are weights [2C 



denotes the projector 



onto the subspace of the total Hilbert space that shall 
never be populated, and gi, minimizes population of this 
subspace (| . The total functional to be minimized, J, is 
given by the sum of , ga and gb, 



D 



9a+gb- 



(4) 



Solving Eq. Q with an iterative procedure in which the 
reference field, £ref{t), is chosen at each level of itera- 
tion to be the optimized field from the previous itera- 
tion ensures vanishing of ga at the optimum, since this is 
reached when the optimal value of J is determined only 
by and possibly gb, but not by the pulse fluence A 
monotonically convergent algorithm is obtained for this 
control problem using a simplified version of Krotov's 
method [1, @ . 

The core idea of Krotov's method 0, Q consists in 
disentangling the interdependence of the states and the 
field. This is achieved by separating the final-time and 
intermediate-time dependencies of the total functional J 
and adding to this a vanishing quantity, ^t—^o~Jq ^dt, 
which is expressed in terms of a functional $ that de- 
pends only on the states and not on the field. The free- 
dom of choice in $ is utilized to ensure monotonic con- 
vergence of the algorithm. Specifically, expanding $ up 
to second order in the states, {ipk{t)}, the expansion co- 
efficients are chosen such that the first and second order 
derivatives of J with respect to the states fulfill the condi- 
tions for the extremum and maximum of J, respectively, 
when J is minimized. Since J thus takes the worst pos- 
sible value with respect to the choice of basis states, any 
change of J due to varying the field, e{t), then leads to 
improvement toward the actual target, minimization of 
J. 

For the simple functional of Eq. ^ , it turns out that 
the second order conditions are trivially fulfilled by cor- 
rect choice of the sign of the weights Xa and Afc, and the 
second order expansion coefficients of $ can be set to zero 
0. The ensuing algorithm thus coincides with that ob- 
tained from straightforward variation of J with a specific 
discretization of the coupled control equations M. It 
yields pulses that implement the desired gate with high 
fidelity, provided that the dynamics allow for it [ll| . For 
functionals with higher than quadratic dependence on 



the states, however, it is essential to include the second 
order contribution to $ in order to ensure monotonic con- 
vergence. The second order expansion coefficient can be 
estimated either analytically or numerically, as detailed 
in Ref. tl2|]. 



B. Geometric theory of non-local two-qubit 
operations 

The group of all two-qubit gates, SU{A), consists of lo- 
cal operations, SU{2)(E) SU{2), and non-local operations, 
SU{A)\SUi2) SUi2). This is a direct resuh of the ex- 
istence of a Cartan decomposition of the corresponding 
Lie algebra T3]. Any two-qubit operation can be written 
as 



U = kiAkz 



where 



exp 



(5) 



(6) 



and the k„ S SU{2) ® SU[2) are local operators. The 
set SU{^)\SU{2) (g) SU{2) is generated by the maximal 
Abelian subalgebra of the su{4) algebra which is spanned 
by the operators <Tj ® <Tj, j = x,y,z. Since these two- 
qubit operators commute, the operations belonging to 
SU{'i)\SU{2) (S) SU{2) can be represented by three real 
numbers (c^;, Cj,, c^). Due to the periodicity of the com- 
plex exponential, the Cj take their value from a three- 
dimensional cube with edges I = [0,7r]. The op- 
erations from this set create and change entanglement 
between two qubits (with the exception of the identity 
operation 1 and the SWAP gate for which all Cj — 
mod TT and all Cj = 7r/2, respectively). Each local op- 
eration ki/2 G SU{2) (g) SU{2) accounts for additional 6 
parameters, yielding the remaining 12 parameters that 
are required to fully characterize the elements of SU{'i). 

The representation of the non-local content of U in 
terms of the Cj is not unique. Different points in the cube 
I'^ may correspond to the same two-qubit operation up to 
local transformations Q. For example, all eight corners 
of the cube are equivalent to the identity operator up to 
local transformations. This symmetry is characterized 
by the Weyl refiection group. It is generated by permu- 
tations or permutations with sign fiips of two entries in 
(ca;, Cj,, Cz). Symmetry reduction of the cube leads to 
a geometric representation of non-local two-qubit gates 
within the Weyl chamber a"*" which is the tetrahedral 
segment of the cube I'^ spanned by {cx,Cy, Cz) = (0, 0, 0), 
(tt, 0, 0), (7r/2, 7r/2, 0), and (7r/2, 7r/2, 7r/2). Ah two-qubit 
gates that are equivalent up to local operations k„ are 
geometrically represented by a single point in the Weyl 
chamber (except on its base where local equivalence 
classes may be represented by two symmetry-equivalent 
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FIG. 1: (Color online) Optimization of non-local gates in 
the Weyl chamber. Local invariants of the B gate (orange) 
and CNOT (red) are approached iteratively (each violet /blue 
point corresponds to one step in the optimization of the ef- 
fective spin-spin Hamiltonian, cf. Section Till Al below) . 

points) Q. For example, CNOT and the controlled tt- 
phase gate (CPHASE) are in the same local equivalence 
class, which is defined by the point (7r/2,0,0), cf. Fig. [T] 
(left panel). 

Each class of all the two-qubit gates that are equivalent 
up to local operations, [U], can also be characterized by 
three real numbers [14| . the local invariants gi, (72 and 

.93- 

It was shown in Ref. [3j that there exists a one-to-one 
correspondence between the points in the Weyl chamber 
{cx^CyyCz) and the local invariants. For a given unitary 
U, 51, (72, 93 are easily calculated from the representation 
of in the Bell basis 0]. Denoting the transformation 
from the logical basis into the Bell basis by Q, Ob = 
Q^UQ = 61F62 where 6„ = Q''^k„Q S SO{i) and thus 
6^6„ = 6„6^ = 1. The matrix F is diagonal. 

Introducing a matrix mu, 

tfi[/ = UbUb = 62F 62 , (8) 

leads to automatic elimination of the first local factor 61. 
The local invariants are defined Q as 

51 = i-«neTr{tfia}' (9a) 
16 

= cos^ Cx cos^ Cy cos^ Cz — sin Cx sin Cy sin Cz , 

92 = ^3mTr{muy (9b) 
Id 

= — sin 2cx sin 2cy sin 2cz , 
g, = ^TriruuV - Tr{rhl} (9c) 

2 2 2 2 2 2 

= 4 cos Cx cos Cy COS Cz — 4: sin Cx sin Cy sin Cz 
— cos 2cx COS 2cy cos 2cz , 













32 


53 


r-n 1 











1 





3 


[CiN U i J 


TV/ A 


U 


U 


U 


U 


1 
i 


[CPHASE] 


Tv/2 














1 


[B-gate] 


Tv/2 


7r/4 














[VSWAP] 


7r/4 


7r/4 


7r/4 





1/4 





[SWAP] 


Tv/2 


Tv/2 


7r/2 


-1 





-3 



TABLE I: Examples of local equivalence classes, their coordi- 
nates (cx, Cy,Cz) in the Weyl chamber, and their local invari- 
ants gi, 52, 93 



where the remaining local factor 62 is eliminated due to 
the cyclic permutation invariance of the trace. To gen- 
eralize the local invariants to the elements of the group 
U{A), i.e., = e*"U' e f/(4) where O' S SU{4), the 
global phase e'" is eliminated by dividing the local in- 
variants by det{U} = e'**". The final form of the local 
invariants is then 

91 = fHeTr{tfi[/}Vl6det{U}, (10a) 

92 = 3mTr{m[/}Vl6det{U}, (10b) 
53 = Tr{my}2- Tr{m2y}/4det{0}. (10c) 

A few examples of local equivalence classes, their coor- 
dinates Ci in the Weyl chamber and the corresponding 
Makhlin invariants gi are given in Table HI 

C. Optimization functional based on the local 
invariants 

To act as a suitable optimization functional, any real- 
valued functional J[{'Pfc}] should fulfill two necessary 
conditions: (i) J must strictly take its global optimum 
for all sets of states {(fik} that represent the desired out- 
come, and (ii) J must be regular, i.e., at least twice 
differentiable. Moreover, as a matter of practicality, it 
should be possible to express J[{95fe}] explicitly in terms 
of the states, {fk}, in order to carry out the differenti- 
ation. While the Weyl chamber coordinates, {cx,Cy,Cz), 
arc only implicitly given in terms of the evolution U of 
the system, the local invariants gi, g2, 93 can directly 
be calculated from U by use of Eqs. ^ (or Eqs. pU]) 
for U G U{4:)), and hence from the time-evolved basis 
states Ifkit)) — 0\(pk{0)). Therefore the local invariants 
lend themselves naturally to the definition of a distance 
measure between the desired local equivalence class [O] 
and the actually realized local equivalence class [U] . We 
define this distance measure as 

3 

d = ^Ag2 with Ag, = |5,(6)-5,(Ut,w)|. (11) 

The distance d constitutes one component of the desired 
optimization functional. Since the time evolution of the 



4 



physical system typically occurs in a Hilbert space that 
is much larger than the logical space, U(r, 0; e) needs to 
be projected into the logical subspace, 

UT,7V = PArU(T,0;e)PAr. (12) 

Correspondingly, d needs to be augmented by a term 
that enforces unitarity in the logical space at the final 
time T, leading to the following optimization functional 
for a non-local equivalence class: 

J|;^ = A.g2 + A.g| + A.g2 + l-lTr{0T,jvUi^} . (13) 

Using the definition of the local invariants, the Agi are 
expressed in terms of the time-evolved basis states as 
follows. Expanding U in the orthonornial basis \ipk), the 



matrix elements of rhu are second order in the states, cf. 
Eq. ([5]). Therefore the gj are fourth order in the states, 
cf. Eq. (|9]), and J^^ turns out to be a polynomial of 
eighth order in the states. The specific form of J^^ is 
given in Appendix [X] 

For such a non-convex functional, the commonly used 
optimization algorithms including that of Ref. [2|] are not 
sufficient to ensure monotonic convergence and indeed 
fail to converge. To the best of our knowledge, Krotov's 
method is the only approach providing an optimization 
algorithm that ensures monotonic convergence for arbi- 
trary functionals in quantum control jl2| . For difficult 
control problems, monotonic convergence is essential to 
reach any optimum, even a local one. Applying Kro- 
tov's method to the functional J^^ , the optimal field is 
obtained iteratively according to 



N 



>fe=l 



de 



(i+i; 



{t) 



Mt) 



de 



,(14) 



where|A</7^'+'^(i)) 



(fi'i^t)). The details of 



the iterative algorithm are presented in Appendix |B] In 
brief, Eq. ([T4)) implies propagating the basis states \ipk) 
forward in time and the adjoint states \xk) backward in 
time with the 'initial' condition, \xk{T)), determined by 
the functional J^^ . S{t)/Xa and a{t) are parameters of 
the optimization algorithm that are constructed following 
the prescription of Ref. [l2| , cf. Appendix IB] 



^1 ^1 ^1/ ^/i ^ 

and combining kj^, kj, kj^, kj yields ki, k2. Assuming the 

errors associated with the local operations ki, k2 to be 
small compared to the error of the non-local operation, 
the actual gate error, £, is then obtained by evaluating 

for kiOT,Ark2, 



£^l-l9.e 



Tr{6^kiOT,wk2} 



D. Actual gate operation and gate error 

Optimizing the functional jji^ yields some gate Ot,n 
that is, up to some small error, in the local equivalence 
class of the desired gate O. In order to actually im- 
plement O, we need to determine the local operations 
ki and k2 such that O = kiUk2- This is achieved by 
transforming both O and Ut,jv to the canonical form A, 

k]^6k2 — A, ki(jT,N^2 = ^- Since A is diagonal in the 
^/ ^1 III ^11 

Bell basis, k^, kg, kj^, k2 are obtained by diagonalization, 
and their combination yields ki, k2. We first determine 
the local operations kj^, kg that transform Ut,jv into the 
canonical form A. This is achieved by calculating the 
gi{^T,N) which yield the Ci and thus A. In the Bell ba- 
sis, A is diagonal and V.^ and k2 are the transformations 

that diagonalize \^T,N^T n ^^-^ Uj. ^UT,Af, each yielding 
- 2 

A . Care must be taken to assure the same ordering of 
eigenvectors when determining V.^ and 1*2. Repeating the 
same procedure for the local operations k^ , k2 that trans- 
form the target operation O into the canonical form A, 



III. APPLICATIONS 

We apply the local invariants optimization functional 
to two examples, an effective spin-spin model that can be 
realized by trapped polar molecules and a Rydberg gate 
for trapped atoms. In the first example, the Hamiltonian 
may become complex, making it impossible to determine 
a priori which two-qubit gates it can implement. The 
Hamiltonian of the second example can realize diagonal 
two-qubit gates only. The complexity that necessitates 
use of optimal control theory in this case draws from 
coupling the logical basis to external degrees of freedom. 

A. Two-qubit gates for an effective spin-spin model 

Trapped polar molecules with ^Si/2 electronic ground 
states, subject to near-resonant microwave driving induc- 
ing strong dipole-dipole coupling, give rise to the effective 
Hamiltonian 

He//(0 = Ml ^ a,A,,{x,,t)a, (15) 



within second-order perturbation theory in the field Q 
(0-4 — 1). The couphngs Ay(a;o,t) = \^{€)\aij{xQ) de- 
pend on the distance between the molecules and on 
the polarization, detuning and possibly time-dependent 
envelope of the microwave field. We consider here SrF 
molecules in an optical lattice with a lattice spacing of 
300 nm and microwave radiation of about 15 GHz. The 
qubit is represented by the spin of the valence electron of 
the molecule in its rotational ground state, as described 
in Ref. Q- We seek to implement two-qubit gates that 
are locally equivalent to the B gate [lH and to CNOT, 
cf. Fig. [1] and Table HI The B gate allows for generating 
a generic two-qubit operation from just two successive 
applications. Arbitrary two-qubit operations can thus 
be implemented with a minimal count of two-qubit and 
single-qubit gates [isj . 

For a single microwave field, it is straightforward to use 
the methods of Section lll DI to determine which non-local 
equivalence classes are accessible under time evolution 
with However, multiple fields are employed, as pro- 
posed in Ref. Q , to allow generation of a broad range of 
effective spin-spin Hamiltonians. Whenever the spin-spin 
interactions deriving from different fields do not com- 
mute, it becomes a non-trivial task to determine which 
two-qubit gates may be efficiently generated by time evo- 
lution under the combined effective Hamiltonians. Op- 
timization of the non-local content of the quantum gate 
reached by time evolution then provides a useful route to 
find the acccessible two-qubit gates, an important task 
for quantum simulations with these effective spin-spin 
Hamiltonians. We illustrate this with an example of de- 
termining which gates are accessible under irradiation by 
two microwave fields with different polarizations. When 
one field is pulsed and the other in continuous wave (cw) 
modality, the effective Hamiltonian is of the form 



H(t) = Ho + ^(i)Hi 



(16) 



where Hq and Hi do not commute, and S'(i) denotes the 
envelope (shape) of the pulsed field, < Sify < 1. We 
choose the polarizations to be ag = 1/a/2, a+ = 1/V2, 
and ~ for the cw field and — 0, a]. — I/V2 
for the pulsed field. We will first optimize for the CNOT 
gate. For this case we assume a rotational transition 
detuning of 1.2 kHz and a Rabi frequency of 590 kHz for 
the cw field and a pulse detuning of 50 kHz. Then the 
drift Hamiltonian, in MHz, in the logical basis becomes 



Hr 



/5.711 
0.324 
0.324 

V 



0.324 
-1.840 
1.054 




0.324 
1.054 
1.840 





and the control Hamiltonian, in MHz, is given by 



Hi 



S{t) 



/-153.65 



V 3.906 





153.65 
16.085 






16.085 
153.65 
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FIG. 2: (Color online) Gate error vs iteration step for direct 
(black dashed lines) and local invariant (red dashed-dotted 
lines) optimization of the CNOT and B gates with time evo- 
lution generated by Eq. (|16|l . The blue dotted line shows the 
gate error for direct optimization of a specific instance of the 
local equivalence class [CNOT] (see text). Red solid curves 
display the value of the local invariants functional J^^ (for 
the direct optimizations the gate error is equal to the value 
of the functional). 



The peak Rabi frequency of the pulse is 1.81 MHz. We 
will then optimize for the B gate. In this case the po- 
larizations and therefore the structure of the Hamilto- 
nian matrix are the same but the numerical values are 
slightly changed since we take different field parameters. 
Specifically, for the cw field we take the rotational tran- 
sition detuning and Rabi frequency to be 1.2 kHz and 
4.74 MHz, respectively, and for the pulse field we take 
detuning and peak Rabi frequency values of 84 kHz and 
1.81 MHz, respectively. Figure [2] shows that direct opti- 
mization for CNOT and B gates is not successful, failing 
to find any high quality solution after a large number 
of iterations (dashed black lines). Thus with the struc- 
ture deriving from this combination of microwave fields 
and polarizations, the Hamiltonian cannot generate the 
unitary transformations corresponding to the B gate and 
CNOT. In fact, it is not evident which gates from which 
equivalence classes can be realized from simply inspecting 
the Hamiltonian. 



In contrast to this lack of success with direct optimiza- 
tion, local invariant optimization of J = J^^ + ga can 
be successfully used to find the time-dependent envelope 
S{t) such that the microwave fields irnplement gates that 
are locally equivalent to the B gate [15] and to CNOT. 
(17) The errors for the optimized gates are shown as a 
function of iteration number in Fig.[2]and compared with 
the results of direct optimization. Progress of optimiza- 
tion in the Weyl chamber is illustrated also in Fig. [T] The 
errors for the Jj^^ optimized gates in Fig. [5] are seen to 
be of order 10~^ (red dashed-dotted line). The result of 
, . standard direct optimization is in stark contrast 

^ ' to this, essentially failing to find a solution, as evidenced 
by the gate error remaining of order unity after 200 iter- 
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ations (dashed black lines). This dramatic difference re- 
flects the fact that the Hamiltonian i^effj when realized 
for the field combination of Eq. (fT6)) . cannot generate an 
evolution that directly yields either CNOT or B gates. 
However, the successful local invariant optimization can 
subsequently be used as input for a direct optimization 
of CNOT, as follows. Inspecting the solution for the local 
equivalence class [CNOT] that is obtained from optimiz- 
ing J^^ shows that in this case the optimal unitary trans- 
formation is a sum of Ud = — diag(l— i, 

which is also in [CNOT] and a smaller Uod with all matrix 
elements except [/23, U32 zero. This provides motivation 
to specify Ud as a target for direct optimization using 
. We find that this leads to a solution of similar qual- 
ity as using J^^ (dotted blue versus dot-dashed red line in 
Fig. In contrast, the optimization of the local equiva- 
lence class [B] does not appear to result in an analogous 
dominant unitary within its equivalence class. Thus it is 
not possible to "guess" which representative of [B] might 
be implemented directly. This example illustrates how 



optimization of the locally invariant functional can be 
used to determine what gate operations are achievable, 
given a possibly intricate Hamiltonian. 



B. Rydberg gate with trapped neutral atoms 

An application to time evolution occurring in a Hilbert 
space much larger than that of the quantum register is 
given by qubits encoded in ^^Rb atoms trapped by op- 
tical tweezers [l6l ]. A non-local gate is implemented by 
simultaneous excitation to a Rydberg state using a near- 
resonant two-photon transition [5| . 

In the experiment of Ref. |l6{ . the qubit states are 
taken to be |0) = |5si/2,F = 2,Mf = 2), |1) = 
\^sif2,F — 1,Mf — 1), the Rydberg state \r) ~ 
|58c?3/2,-F = 3,Mp — 3), and the intermediate state for 
the two-photon transition |i) — \5pi/2,F — 2,Mp = 2). 
In the rotating- wave approximation, the Hamiltonian for 
a single trapped atom reads 



(t + y*™p(x,)) +|1)(1|® (t + y,™p(x,)) +|z)(z|® ft+^j +|r)(r|® ft+^ 



(|0)(z| + |z)(0|) 1,^. + ^ (|z)(r| + \r){^\) ® 1,^ 



(19) 



Here, Xj denotes the position operator of atom j, T the 
kinetic energy operator, and rii{t) the time-dependent 
Rabi frequencies of the red and blue lasers {ujr = 795 nm 
and lob = 474 nm, respectively). The maximum Rabi 
frequencies are taken to be fli,o = 2tt ■ 260 MHz, i.e., 
flR_o is equal to and ils.o larger by a factor of 10 than 
those of [lal- ^iit) is parametrized as 

n,{t) = 17,,o(tanhe,(i) + l)/2 e [0, f],,o] , 

with Si{t) determined by optimal control. The detuning 
Sr of the red laser is chosen to be Sr — 2tt ■ 600 MHz, 
slightly larger than in jl6| . The two-photon detuning 
from the Rydberg level is given in terms of the Stark 

shift, 6b = {nl f, - nl ,,)/i4SB,) = 0. 

The total two-atom Hamiltonian includes the long- 
range interaction when both atoms are in the Rydberg 
state, 

H^^^i) = Hi^\t) (g) 14,2 «) +14,1 «)H2^^(t)«) 1x1 

+ \rr){rr\(S)^, (20) 
f 

with r — |xi — X2I the inter-atomic distance, uq is chosen 
to reproduce the estimated interaction energy of 50 MHz 
at ro = 4/zm (mq = 3.284 • lO^E^a.^). With this interac- 
tion, the atoms need to spend a minimum of 10 ns in \rr) 



to pick up a non-local phase of tt. In all other internal 
states, the interaction between the atoms at a separation 
of 4 fim is negligible. We approximate the optical tweez- 
ers trap by harmonic potentials and integrate over the 
center-of-mass motion. The trap of width, a = 0.75 /zm, 
and depth, T^nin = — 4.5/cBmK, is slightly stronger than 
in Ref. fl6\. 

' (2) 

The Hamiltonian, H (t) , is represented on an equidis- 
tant Fourier grid extending for ±0.3 /im around tq. In 
order to evaluate Eq. ([H)) . the Schrodinger equation is 
solved simultaneously for initial states ■ (po(r) — 
0, 1) with ipo(r) the ground state of the trap, using a 
Chebychev propagator. 

The errors from the optimal gates are shown in Fig- 
ure [3] as a function of the corresponding optimal gate du- 
ration T and illustrate several key results. First, the large 
error resulting from direct optimization of CNOT (blue 
- (2) 

triangles) shows that H (t) cannot generate CNOT di- 
rectly. Second, the minimum errors of the Rydberg gate 
are comparable for optimization with and J^^ (with 
a slight advantage for the functional based on the local 
invariants [2lj ) and both of these values reflect the quan- 
tum speed limit [l7| . Third, if spontaneous emission from 
\i) is neglected (filled symbols), the minimum gate dura- 
tion of ~ 40 ns to achieve high fidelities is determined pri- 
marily by the interaction strength in the Rydberg state. 




gate duration T ( ns ) 



FIG. 3: (Color online) Optimal gate error as a function of the 
gate duration for Rydberg gate optimization based on local 
equivalence classes (Jt^) and on direct optimization (Jt) of 
two specific two-qubit gates (yr-phase = diag(— 1, 1, 1, 1) and 
CNOT) without (a) and with (b) an additional constraint 
suppressing population in the intermediate state |i) — |5pi/2)- 
Dashed lines with open symbols (solid lines with filled sym- 
bols): simulation including (neglecting) spontaneous emission 
from Large filled symbols: calculations with trapping po- 
tential kept on during the gate. 

Motion of the atoms in the trap leads to larger minimum 
gate errors for increasing gate duration. Gate errors close 
to the fault tolerance threshold of 10~^ are obtained only 
when the trapping potential is kept on during the gate 
(large filled symbols in Fig. |3ja)). 

The main limiting factor for a high-fidelity implemen- 
tation of the Rydberg gate using this particular near- 
resonant two-photon transition is due to spontaneous 
emission from the intermediate state |z) as is evident 
from Fig. [3] (open symbols). Imposing an additional con- 
straint suppressing population oi\i), cf. Eq. dSj, leads to 
improved solutions but leaves the minimum errors still 
two orders of magnitude above the fault tolerant thresh- 
old (Fig.lSjb)). Near resonant intermediate states should 
therefore be avoided by suitable choice of Rydberg and 
intermediate states, for given laser frequencies. 

Sample Rabi frequencies for optimized pulses and the 
corresponding two-qubit state dynamics are shown in 
Fig. m Without spontaneous emission, the minimum 
gate duration for which a high-fidelity implementation of 
CNOT can be obtained is 40 ns. The optimal pulses con- 
sist of a sequence of fast switches pumping population in 
a ladder-like fashion. This is reflected by the small circles 
of |01) , 1 10) states in the complex plane, indicating clearly 
non-adiabatic dynamics. The time between switches cor- 
responds to the duration of a 7r-pulse at the maximally 
allowed Rabi frequencies fli,o. Correspondingly, the spec- 
tra (not shown) consist of modulated sidebands spaced 
by fii,o- Experimentally, switches of the field amplitude 
on a nanosecond timescale can be generated using stan- 
dard electronics (electro- or acousto-optical modulators). 
When spontaneous emission from the intermediate state 




time ( ns ) Re [ >^'(t) ] 



FIG. 4: (Color online) Left: optimal Rabi frequencies with- 
out (a) and with (b) spontaneous emission from Black 
dotted line represents the initial "guess" pulse form. Right: 
corresponding dynamics of the two-qubit basis states in the 
complex plane. 



is taken into account in the optimization using the addi- 
tional constraint gb, the smallest gate errors achievable 
for the given detunings and maximum Rabi frequencies 
are of a few percent. This requires gate durations of 
50 ns or more, which is larger than the gate durations 
without spontaneous emission. While surprising at first 
glance, it refiects the fact that the only way to avoid 
populating \i) is by adiabatic passage, and adiabaticity 
requires sufficient time. Correspondingly, the optimal 
Rabi frequencies in Fig. Hfb) show a double STIRAP- 
like behavior [31 : In the first half of the time interval, 
the blue pulse connecting \i) and |r) takes the role of 
the Stokes laser and precedes the red pulse connecting 
|0) and \i) which corresponds to the pump laser. In the 
second half of the time interval, these roles are reversed. 
The optimal entangling strategies in the two examples 
are thus based on well-known, robust and feasible control 
schemes - population transfer to and from the Rydberg 
state via 7r-pulses or STIRAP. The additional twist that 
is afforded by optimal control is speed-up, i.e., implemen- 
tation of the shortest possible gate duration. This comes 
at a comparatively low price - additional modulations in 
the optimal fields on a nanosecond timescale translating 
into spectra whose bandwidth never exceeds a few GHz. 
It should therefore be comparatively straightforward to 
implement these optimal pulses in an experiment. 

Both with and without spontaneous emission, the 
qubit phase dynamics (right panel of Fig. 2]) show 
achievement of the desired non-local phase, x = 0oo ~ 
001 "'/'lo+'/'ii I with complete freedom in the local phases. 
This confirms the fact that J^^ imposes fewer constraints 
than direct optimization. Since spontaneous emission 
from intermediate states can be eliminated by suitable 
choice of atomic states and exciting lasers, it is evident 
that this non-local optimal control allows a high-fidelity 
implementation of the Rydberg gate to be achieved even 
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for a setting where the blockade regime is not reached and 
when additional entanglement between qubit and exter- 
nal degrees of freedom is allowed. 



IV. SUMMARY AND CONCLUSIONS 

By constructing a new optimization functional based 
on the local invariants of two-qubit operations, we have 
shown how to extend optimal control theory to take spe- 
cific requirements of quantum information applications 
into account. In particular, we have developed an au- 
tomated way to determine for a given physical system 
whether a desired non-local content can be realized, and 
if so which two-qubit operation in that desired local 
equivalence class can best be realized. We have illus- 
trated the power of this approach with two examples 
relevant to quantum simulations and quantum compu- 
tation. The first example addressed the time depen- 
dent control of effective spin Hamiltonians generated by 
trapped dipolar molecules. It showed that optimization 
of non-local content of two-qubit operations by local in- 
variant optimization allows determination of which gates 
may be reached from complex Hamiltonians with non- 
commutative time dependence. The second example ad- 
dressed the performance of a Rydberg gate with trapped 
atoms. It demonstrated that use of optimal control the- 
ory with the local invariant functional yields a faster gate 
than direct optimization, in the presence of both cou- 
pling between internal and external degrees of freedom 
and spontaneous emission. 

The new optimization functional can easily be 
adapted to target all perfect entanglers, using the ge- 
ometric definition of these derived in Q. A related 
question of optimizing multi-partite entanglement was 
recently addressed using time-local control theory [l9| . 
Local invariant optimization is expected to display its 
full potential in the presence of decoherence, particularly 
when distinct two-qubit gates in the same local equiv- 
alence class are differently affected by the decoherence. 
The required extension to open quantum systems is the 
subject of future work. 
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Appendix A: Explicit form of the local invariants 
functional 



Using the definitions of tfi and the local invariants, 
Eqs. ^ and (fTU|) . the functional is expressed in terms of 
the states, 

Jt = /f + /I + /3 + /I + 1 - ^ Tr { UT,ivO+ ^} 
with 



3det{uT,A.}]-l^ 

k,l 

-2dl0f - 4 (Sfc • /3fc) (a, • a) 
/2 = 3m aodet|uT,Af| - 



-4/3^ [di-Pi 



4afe [di- Pi 



dldf + ^Jf ~ 2dl0^ 



■2 a2 



6odet{UT,Ar} 

-4 {dk ■ h) [m ■ a) - {dk -dif - (/3fe • a)' 
+2 {dk ■ di) (/3fc • a) + 4 {dk ■ di) (pk ■ a) 

U = am [bodet [(Jt,n}] ~ ^Yl ^"fc ' 

k,l L 

-4/3^ (a,-A) -4(afc.«0(afc-A) 

+4 (/3fc • a) (Sfc • a) , 

where the sum runs over the N logical basis states. 
The constants uq, bo are obtained by calculating m, cf. 
Eq. ([5]), for the target gate O, 

V {mo} 



ao 



bo 



16det{0} 

[Tr^{mo}- Trjmg,}] 
4det{6} 



and dk {Pk) is the vector containing all real (imaginary) 
parts of the expansion coefficients of the state \ipk) with 
respect to an orthonormal basis, spanning the 

complete Hilbert space, 

{ak)m = *He [(m|(/3fe(i))] , m = 1, . . . , dim('H) 

{Pk)m = am [(m|v5fe(t))] , m 1, . . . , dim('H) . 

Note that J^-^ is a polynomial of 8th degree in the states. 

Appendix B: Outline of the optimization algorithm 

The algorithm is determined by the functional Jt, ad- 
ditional time-dependent constraints ga, gb and the equa- 
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tions of motion [1^ and given by the following set of 
equations: 

1. Forward propagation to obtain the new states 

|l^rHi))=-^H[e^^+^^]|^^'^W). (Bl) 

The initial states are given by the basis expansion 
of the time evolution operator, i.e., the logical basis 
states in our case. 

2. Backward propagation to obtain the adjoint states 
\Xk\^))' containing an inhomogeneity if gb ^ 0, 

jy:\t))^~l^[^^'^]\x'-:\t)) + \rj) (B2) 

with the 'initial' condition at time t — T deter- 
mined by Jt, 

|Xfc(r))-V(^|JT, (B3) 
and the inhomogeneity 

l^>=V(^|5fc. (B4) 



3. The equation to determine the new field from 
\4^^\t)) and \xf{t)), Eq. (HI, with a{t) given 

by 



a{t) = C{T -t)- A. (B5) 



The constants A and C are parameters of the al- 
gorithm that can be estimated analytically (based 
on a worst case scenario) or numerically |12l |. For 
the local invariants functional, J^^, the analytical 
estimate yields A = 90 for Hamiltonian and 
A = 580 for Hamiltonian (pn|) . while numerically 
A = b and A = 20, respectively, turned out to be 
sufficient. C = for c/b = and C < -Xb/NT for 
9b ^0. 



Compared to a linear version of the Krotov algorithm [2'] , 
the additional effort consists only in storing the forward 
propagated states from the previous iteration, \(p'}!\t)), 
and calculating \Aipk{t)) and a{t). 
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